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Abstract 

We study the renormalization of the properties of low lying charm and hidden 
charm scalar mesons in a nuclear medium, concretely of the -D s o(2317) and the theo- 
retical hidden charm state X(3700). We find that for the D s q(2317), with negligible 
width at zero density, the width becomes about 100 MeV at normal nuclear matter 
density, while in the case of the V(3700) the width becomes as large as 200 MeV. 
We discuss the origin of this new width and trace it to reactions occurring in the 
nucleus, while offering a guideline for future experiments testing these changes. We 
also show how those medium modifications will bring valuable information on the 
nature of the scalar resonances and the mechanisms of the interaction of D mesons 
with nucleons and nuclei. 

1 Introduction 

The modification of the properties of elementary particles in nuclei is a rich field that 
helps simultaneously to learn about excitation mechanisms in the nucleus as well as 
properties of the elementary particles [1J. Important features about the nature of 
certain particles can be better observed in nuclei. Take for instance the A(1520). 
This resonance couples strongly to 7rS(1385), but this is hardly visible in the decay 
of the A(1520) particle since there is no phase space for it. When this resonance is 
placed inside a nucleus the pion can become a particle hole excitation {ph) and one 
gains 140 MeV of phase space for the decay, which results in a considerable increase 
in the width of the A(1520) in the medium [2]. Take another example: the u meson. 
This meson decays into 3tt, which is supposed to go through ptr, but there is no 
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phase space for decay in this channel except for the p width. Once again the tttt 
decay channel of the p will be modified in the nucleus, as well as the individual n, 
which can become again a ph excitation producing a much larger phase space for the 
decay of the u. It is then spectacular to find that the width of the uj in the medium 
at nuclear matter density becomes as large as 100 — 150 MeV OH]. This discussion 
is only qualitative and more subtle details must be considered as done in other works 



Among so many examples in the literature, the case of the renormalization of the 
properties of the scalar mesons has played a special role. One reason is that there is a 
long debate on the nature of the scalar mesons, with different assumptions about their 
nature as qq states, KK molecules, mixtures of qq with meson-meson components, 
or dynamically generated resonances from the interaction of coupled channels of 
two pseudoscalars [8j-[26j. We will study here the properties in the medium of two 
dynamically generated states, meaning states made out of two mesons. Pioneer work 
in the study of dynamically generated resonances was done in [HJ [9] , where starting 
from one seed of qq, a large meson cloud was demanded to explain data of the 
low lying scalars. Further work in this direction followed using the unitary coupled 
channel chiral approach in [241 [25| 126]. In this sense, the predictions of the different 
models on the medium modification of these resonances are important in order to put 
stringent constraints on the different assumptions about the nature of the resonances. 

Among the low lying scalars, the <t(600) has been the most studied. Several 
theoretical approaches have predicted strong medium effects on the tttt interaction 
in the scalar isoscalar (a) channel. In [27], the a and ir mesons at high density 
are studied within the Nambu-Jona-Lasinio model and the authors find that the mo- 
drops subtantially with the density, whereas the m n increase at higher densities. The 
7T was previously studied at finite temperature and density within the same model 
in [28]. In Ref. [29], Hatsuda et al. studied the a propagator in the linear a model 
and found an enhanced and narrow spectral function near the 2tt threshold caused 
by the partial restoration of the chiral symmetry, where m a would approach m v . 
The same conclusions were reached using the nonlinear chiral Lagrangians in Ref. 
[30] . Similar results, with large enhancements in the tttt amplitude around the 2ir 
threshold, have been found in quite different approaches by studying the s— wave, 
1 = tttt correlations in nuclear matter [31] . In these cases the modifications of the a 
channel are induced by the strong p— wave coupling of the pions to the particle-hole 
(ph) and A-hole (Ah) nuclear excitations. A more recent and updated theoretical 
work on the evolution of the a poles in the medium can be seen in [32] . 




On the experimental side, there are also several results showing strong medium 
effects in the a channel at low invariant masses in the A(ir, 2ir) [33] and A(*y, 2tt) [33] 
reactions, which have been the object of study in [35 j . [36] and [37] . 

The /o(980) and ao(980) scalar mesons have also been analyzed in the nuclear 
medium [38], but unlike the case of the a, no experimental action has been taken in 
this case. 

In the present paper we retake research along these lines and study the medium 
modification of the scalar mesons which are dynamically generated in the charm sec- 
tor. Concretely, we shall study the medium modification of the D s o(2317), which is 
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obtained in the theoretical studies of [39j-[42j. Within the context of lattice gauge 
theories, in [43] the authors find hints that there is a DK bound state that can be 
identified with the D s o(2317). In addition, we shall also study the medium modifi- 
cations of a hidden charm scalar meson predicted in [42], for which there are some 
indications that could have been observed in the experiment of the Belle collabora- 
tion [43] through the reanalysis done in [45] . In this experiment a broad bump is seen 
in the DD mass distribution around the DD threshold, which in [45] was shown to 
be better explained by the X(3700) pole of the DD amplitude below threshold than 
by the Breit Wigner distribution proposed in the experimental paper. However, it 
should be noted that a better fit to the data with a Breit Wigner distribution than 
the one of [44] was obtained in [45] . yet with a x 2 value slightly worse than the one 
obtained with the pole below threshold. 

2 Brief discussion on the dynamical generation 

of the L> s0 (2317) and X(3700) 

We follow the details of [42] , where a Lagrangian is taken for the interaction of two 
pseudoscalar mesons. The Lagrangian is an extrapolation to SU(4) of the SU(3) 
chiral Lagrangian used in [241 [26] to generate the scalar mesons <t(600), /o(980), 
ao (980) and k(900) in the light sector, but with the 5f7(4) symmetry strongly bro- 
ken, mostly due to the explicit consideration of the masses of the vector mesons 
exchanged between pseudoscalars in the equivalent theory using the hidden gauge 
formalism for the vector mesons [46] -[49]. A different breaking of SU(A) is also con- 
sidered in [32], following general rules of SU(n) breaking [50], which serves as an 
indication of theoretical uncertainties. We follow here the version based on the hid- 
den gauge formalism, including the exchange of heavy (charmed) vector mesons in 
the Lagrangian. 

We would like to put the <Si7(4) breaking in a certain context. The basic assump- 
tion underlying [42] 143] is that the vertices in the hidden gauge formalism, of four 
vectors or three vectors, are approximately SU(A) symmetric. As a first step, the 
kernel of the Bethe Salpeter equation (the potential) already breaks SU(4) symme- 
try in the terms that exchange heavy vectors, as we have mentioned. Later on, the 
amplitudes calculated with the unitary approach break 5C/(4) symmetry because the 
physical masses of the particles are used to respect exactly thresholds and unitarity 
in coupled channels. This situation is already present in SU(3). The starting lowest- 
order chiral Lagrangians are SU(3) symmetric and the amplitudes obtained break 
SU(3) symmetry due to the different masses of the particles belonging to the same 
SU(3) multiplet. As an example, one of the A(1405) states [51] and the A(1670), 
which appear in the approach of [51, 52J, are degenerate if the masses within the 
same SU(3) multipletes are taken equal. One can see there an example of a large 
SU(3) breaking present in nature, which, however, is consistent with assuming an 
exact SU(3) meson-baryon interaction Lagrangian. 

With these assumptions for the SU(4) symmetry and its breaking one gets realistic 
results for the spectra of mesons in [42] 145] , which also agree with those obtained in 
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the heavy quark formalism for the case of light-heavy meson interaction [53], up to 
a mass term with no practical consequences. It is also interesting to note that the 
same basic SU(4) assumptions are done for the interaction of mesons with baryons 
in [54l ESI |56] and one obtains realistic results concerning the A c (2595) and £ c (2800) 
resonances. 

Following [24] , one derives the kernel (potential) from the lowest order Lagrangian 
and iterates it to generate all the terms of the Bethe Salpeter series, which can be 
summed up in the on shell formalism |57} [58] by means of 

T=[1-VG]~ 1 V, (1) 

where T is a matrix in the space of coupled channels representing the transition 
scattering amplitude from one channel to another and V the equivalent matrix for 
the transition potential. The diagonal matrix in the coupled channel space G accounts 
for the loop integral of the two particle propagator of any intermediate state 
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where P is the total four momentum, q one of the pseudoscalar four momentum, p 
the on shell three momentum and mi, 777,2 the masses of the two pseudoscalar s. 

This integral requires dimensional regularization by means of a substraction con- 
stant, aj, or cutting the three dimensional integral in G with a cut off. Both methods 
establish equivalent schemes in a certain chosen region of energies [58] . However, the 
use of the dimensional regularization method relies upon Lorentz covariance of mag- 
nitudes, which is lost in nuclei where one has a privileged reference frame, the one 
where the nucleus is at rest. As a result, the use of the dimensional regularization 
method introduces pathologies when including the selfenergy of the particles in the 
medium, which are avoided with the use of a cut off [59J. Hence, in the channels 
where we renormalize the particles in the medium, we stick to the cut off formalism. 
For all the other channels we use the same dimensional regularization approach of 
|42j which guarantees unitarity. The cut off method is used only in the channels DK 
and DD, and the values of the cut off will be shown later on, but they are much 
bigger than the on-shell three momenta of the particles in the loops, such that the 
imaginary part of Gu is obtained exactly in the free case. In free space we have 

q ( s \ = f qmax 1 2d 1 ^1 +^2 (4) 

where q ma x stands for the cut off, Ui = (qf + m 2 ) 1 / 2 and the center-of-mass energy 
(P°) = s. The expression of Ga will be changed to account for the medium effects 
on the pseudoscalar mesons in Sec. HI 
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The T matrix of Eq. ([I]) generates poles for some quantum numbers. We look for 
them in the second Riemann sheet for the channels which are open and those poles 
are associated to resonances. For the channels where we found bound states, i.e., 
states below threshold, the poles appear in the first Riemann sheet. This is actually 
the case for X(3700) in the DD and also for the Dq s {22>\7) in the DK channel. Close 
to a pole, the amplitude looks like 



where Re zr gives the mass of the resonance and Im zr the half width. The constants 
gi, obtained from the residues of the amplitudes, provide the coupling of the resonance 
to a particular channel and indicate the relevance of this channel in building up the 
resonance. This said, it is useful to recall that the X(3700) and D s o(2317) are 
obtained from poles of the scattering matrix with quantum numbers (C = 0, S = 0, 
1 = 0) and (C = 1, S = 1, I = 0), respectively. In Tables 1 and 2, we show 
the coupling of each state to the different channels contributing to these quantum 
numbers, which are obtained from [32] but including the if and the mixing with i] 
(see [60]). 



Channel 


Re(g x ) [MeV] 


lm{g x ) [MeV] 


\9x\ [MeV] 


7T + 7T~ 


9 


83 


84 


K + K~ 


5 


22 


22 


D+D- 


5962 


1695 


6198 


n 7i° 


6 


83 


84 


K°K° 


5 


22 


22 


f]rj 


1023 


242 


1051 


r\rf 


1680 


368 


1720 


rfrf 


922 


-417 


1012 


D°D° 


5962 


1695 


6198 


DtD- 


5901 


-869 


5965 


VcV 


518 


659 


838 


VcV' 


405 


9 


405 



Table 1: V(3700): Couplings of the pole at (3722-il8) MeV to the channels (C=0, 
S=0,I=0). 



From these tables we observe the following: 

1) The heavy singlet, hidden charm state V(3700) couples most strongly to DD. 
Next it couples to D S D S . However, from the square of couplings which would enter 
into the selfenergy of the resonance due to the intermediate meson-meson states, one 
gets a factor of two more weight for the DD than the D S D S states. On the other 
hand, while the pole around 3700 MeV is close to the DD threshold of 3738 MeV, 
it is about 238 MeV away from the D S D S threshold of 3938 MeV. The off shellness 
of the D S D S in the loop function of Eq. further reduces the contribution of this 
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Channel 


Re(g Ds0 ) [MeV] 


lm(g Ds0 ) [MeV] 


\9dJ [MeV] 


K+D° 


5102 





5102 


K°D + 


5102 





5102 


V D+ 


-2952 





2952 


rfD+ 


4110 





4110 


VcDt 


2057 





2057 



Table 2: D s0 (2317): Couplings of the pole at 2317 MeV to the channels (C=l, S=1,I 



channel to less than 10% of the DD. This is important to note since, when dealing 
with the medium correction in the next sections, we shall consider the normalization 
of D but not the one of D s . 

2) The -D s o(2317) couples most strongly to DK. The next channels are the D s r/' 
and D s r/. The same considerations as before lead to a similar relative contribution 
of D s r]' and D s rj with respect to the dominant DK channel. Once again, we are 
lead to deal with the medium modification of the D meson, those of the K state 
being relatively unimportant since the KN interaction is not too strong and has no 
singularities [61]. 

It should be noted that in the step from SU(3) to 5C/(4) we introduce more 
uncertainties than one has in SU(3). That is the reason why in [42] two different 
models, which break SU(4) in different ways, were used to see the uncertainties 
in the results. It was found there that the results for the V(3700) and _D s o(2317) 
states were rather independent of the model, while other states predicted here were 
more model dependent. Even though, it is important to make an evaluation of the 
theoretical uncertainties for those states and, for this purpose, we have followed the 
same approach as in [32] • Thus, we evaluate the results using a sample of parameters 
which are varied within the reasonable limits discussed in [42 . In particular we vary 
the fn and fo parameters in the range /„- £ [85, 115] MeV and fr> 6 [146, 218] MeV. 
This alone gives a good approximation to the uncertainties and we will perfom the 
medium calculations evaluating the uncertainties with this method. 

The main results of the paper are the effects of the medium in the D s o(2317). 
The results for this resonance were obtained in [32] , where the prescription of using 
/tt for the light mesons and fjj for the heavy ones was used. For the case of scattering 
of light mesons with heavy ones, chiral symmetry requires the use of f n in all cases 
(eventually fx if kaons are involved). For this reason we redo the calculations with 
the new chiral symmetric prescription. We also introduce the novelty with respect 
to [32] of considering also the rf in the set of pseudoscalar mesons according to the 
method of [60]. We should bear in mind that the substraction constant of Eq. (|3|) 
was slightly tuned from its natural value to get the mass of the Z) s o(2317) at the 
right place. Hence, we do the same here and then we look at the results obtained for 
the couplings to the channels in the two cases. Those are shown in Table El As we 
can see, the results are very similar. For the most important building block, the DK 
channel, the differences of the coupling are of 4%, indicating that the errors induced 
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Channel 


\gi\ Model A [MeV] 


\ gi \ Model B [MeV] 


DK 


7215 


7503 


D s r] 


2952 


3005 


D s T]f 


4110 


4146 


DsVc 


2058 


1246 



Table 3: Couplings of the _D s o(2317) to its building blocks. Model A refers to the model 
using both, /„. and fjj in the couplings, while in Model B only /„. is used, respecting 
constrains from chiral symmetry. The channels are in isospin basis. The position of the 
pole is fixed in both models to 2317 MeV, taking an = —1.48 in the model A, and 
an = —1.16 in the model B (a# means the substraction constant used in [12] for the 
channels involving at least one heavy pseudoscalar meson). 



by the explicit chiral symmetry breaking of [42] are very small. Yet, in the present 
paper we shall use the chiral symmetric version described here. Using fx instead of 
/tt leads to much smaller differences than those in Table [3l However, we shall use 
just f n . 

For the X(3700), which comes mostly from DD, we have no such constraints 
from chiral symmetry and we follow the approach of [32], except for the inclusion 
of the rf channel. We evaluate uncertainties in the results, though, by using again 
/tt instead of both fr> and /„-. The stability of the couplings with respect to the 
changes done can be partly justified by recalling that, with one channel dominance, 
the coupling of a bound state to its constituents is given in terms of the binding 
energy by the compositness condition of Weinberg [62]-[67]. In the present case, the 
DK is the dominant channel, but other channels also matter. This is why some 
changes, although small, were found in the couplings while demanding the same 
binding energy. 

3 The selfenergy of the D meson 

3.1 s-wave selfenergy 

We shall use the T = results from the work of [56]. There, the D meson selfenergy 
is obtained from a selfconsistent coupled channel calculation, whose driving term is a 
broken SU(4) s-wave Weinberg Tomozawa interaction supplemented by an attractive 
isoscalar-scalar term. The introduction of a supplementary scalar-isoscalar interac- 
tion in the diagonal DN channel, the S_djv term, which is prevalent in the QCD 
sum rule and mean-field approaches, is, however, a subject of controversy. Its effect 
on the D meson self-energy was studied in [55] and [56], and compared to the case 
where this term was neglected. It was found in [56] that the results obtained with 
or without this term of uncertain origin were qualitatively identical, such that the 
phenomenology did not allow to draw any conclusion concerning it. The differences 
found there between the two options are far smaller than the uncertainties that we 
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have in the present problem from other sources and, hence, we take the option of 
ignoring this term. The Bethe Salpeter equation is then solved using a cut off reg- 
ularization, which is fixed by reproducing the position and width of the A c (2593) 
resonance. As a result, a new resonance in the 1 = 1 channel is generated around 
2800 MeV, the S c (2800). The coupled channel structure includes: vrA c , vrS c , DN, 
rjA c , KE C , rjT, c , KE' C , D S A, D S T,, r/A c and rj'E c . 

The in medium s-wave DN amplitude accounts for Pauli blocking effects on 
the nucleons in the DN channel, mean-field bindings of baryons via a a-u model, 
and renormalization of it and the D through their corresponding selfenergies in the 
intermediate propagators. The s-wave D selfenergy is obtained iteratively following 
a selfconsistent procedure as one integrates the in medium s-wave DN amplitude 
over the nucleon Fermi sea n(p): 

U ( g(q°,q,p) = J ^n(p)[f^ \p ,P,p)+3f ( J; 1 \p°,P J p)} , (6) 

where T^ ' 1 ) stands for the in medium s-wave DN amplitudes in / = and 1=1. 
The quantities P° = q° + E^{p) and P = q+p are the total energy and momentum of 
the DN pair in the nuclear matter rest frame, with E^{p) being the single-particle 
nucleon energy and the values (q°, q) the energy and momentum of the D meson 
also in this frame. For more details see [56j . 

In fact, the s-wave D meson selfenergy in the medium was initially studied in 
|68[ [69l [M] and further work was done in |55|, [56] . There are novelties in the approach 
of |55[ [56] with respect to the earlier works. Indeed, for consistency with the reduc- 
tion from t-channel vector-meson exchanges to a zero-range Weinberg- Tomozawa 
form, corroborated with explicit cancellations of terms, the model of [551 156| removes 
a factor k^k^/My, which was used in [53], and by means of which a better width for 
the A c (2593) is obtained in [551 [56] . Also in [551 [56] the authors use a conventional 
momentum cut-off regularization that was found to be more appropiate than the 
dimensional method in view of its application to meson-baryon scattering in the nu- 
clear medium where Lorentz covariance is manifestly broken, as discussed in Section 
2. 

With respect to the work of |68[ [69] there are also novelties in |55[ 156] . In the 
exploratory work of [681 [69], the free amplitudes were constructed from separable 
coupled channel interactions obtained from chiral motivated Lagrangians upon re- 
placing the s quark by the c quark. While these works give the first indication that 
the A c (2593) could have a dynamical origin, they ignored the strangeness degree of 
freedom due its very construction. Therefore, the ir and K (Goldstone) mesons were 
not treated on an equal footing, and the role of some channels that would appear in 
the corresponding SU(4) meson and baryon multiplets was ignored. 

3.2 j9-wave selfenergy 

We start by recalling the SU(3) chiral Lagrangian [70l [71] for the coupling of pseu- 
doscalar mesons of the octet of the tt to the baryon octet of the proton p 

4 B) =\d< 57 m 75K, B) > +h < B^K, B] > , (7) 



S 



where 



u 2 = U = e l — , (8) 



with <p the usual SU(3) matrix of the meson fields, / = 1.15/^ with f n = 93 MeV 
and 

= ivtdpUv! = -1- drf + O(0 3 ) . (9) 

The B and -B terms stand for the SU(3) matrices of the baryon fields and <> 
for the trace in SU(3). Hence, at the one meson field level we have 

£ i B)= -7Tf D <B^{d^B}>-^- f F< B^ ls [d^,B] > . (10) 

The u{p')^ tJ "y^u{p) vertex, assuming p ~ since it will be the momentum of a 
nucleon in the Fermi sea, can be expressed up to 0(1/M 2 ) in terms of the a operator 
such that the T matrix corresponding to the diagram of Fig. [1] is given by 

\ 

X 

\ 

, \ , 



p p' 



Figure 1: Meson-baryon scattering with an outgoing baryon. The labels p, p' and q refer to the 
momenta of the initial baryon, final baryon and meson respectively. 



it 



V2f 



a-q[l 



2M' 



[(D + F)<Bcf)B>+(D-F)<BB(t)>} , (11) 



with M' the mass of the outgoing baryon in Fig. [TJ We take D = 0.80 and F = 0.46 
from [72l E3J E3] • in order to evaluate the coupling of the D meson to the nucleon 
and A c , S c we use SU{4) symmetry. We couple the 20-plet of the baryons, to which 
the nucleon belongs, to the 20 representation of the antibaryons in order to give the 
15-plet of the mesons of the ir and the D [75]. By using the SU(4) Clebsch-Gordan 
coefficients of [76] , we have two independent irreducible matrix elements which can 
be related to the D and F coefficients. The result is that the couplings D p — > A+, 



D°p -> 
of K~p 
[77] by 



-, D+p -> 
A, K~p 



D°n - 
, K°p 



y0 



D + n - 
, K~n 



A 



D + n — > are identical to those 
", K°n — > A, K°n — > S° given in 



iVdny = ^ 



1 



2M' 



D + F D 
a — — h p 



F 



V 



V 



(12) 



with the coefficients a, (3 of the Table H] 

We also take into account the coupling of the D meson with S*(2520) and N, in 
analogy to the p-w&ve interaction of pions and kaons with nucleons. For pions and 
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D°p -> A+ 


D°p -> E+ 




D+p -> £++ 


a 
P 


•2 

V3 
1 

^/3 



1 


"73 
^ A "I 









Table 4: 


Coefficients for the DNY couplings 





kaons it was shown that the iV X A and A r 1 S*(1385) excitations, respectively, are 
relevant for the calculation of the p-wave self-energy. Once again, we obtain the same 
result as in [77J for the iV- 1 E*(1385) 

-Wdny* =aS^ -q ^~^~~^~J ' ( 13 ) 

with being the spin 1/2 — ► 3/2 transition operator and a the coefficients given in 
Table El 





D°p -> £*+ D°n -> £*° D+p -> £*++ D+n -> E*+ 


a 


_ i _i _! i 



Table 5: Coefficient for the DATE* (2520) couplings 




Figure 2: p-wave selfenergy diagram of the D meson. 



With all those couplings, we can readily evaluate the p-wave D selfenergy given 
by the diagram of the Fig. [21 in complete analogy to [77] . The p-wave contribution 
coming from the N~ 1 A C and N Yl c excitations reads 

+ \Bl 0n ^q 2 U n {q\q,p)}Fl{q 2 ) , (14) 

where 

Bdny = ( 1 - m^) 



D + F D 
Ot—r-„ V P 



2/ 



2/ 



(15) 
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and U is the Lindhard function for the iV X Y excitation given by 

/ n s 3 My f 1 , 9n I 2 + 1 I! 

Im = -""TP— {(! - z 2 )0(l- I 2 |)} (16) 

My 




(My - M) 



QPf 



with p = p n + p p , the nuclear density, f>p = (3/T 2 p/2) 1 / 3 the Fermi momentum, My 
the hyperon mass and M the nucleon mass. The same result holds for the p-wave 
D + selfenergy ignoring small mass differences between particles of the same isospin 
multiplet. 

The p-wave selfenergy due to the excitation of the decuplet is also readily evalu- 
ated and we find 

ngV.O) = {\c 2 D0p ^q 2 U^(q\q,p) + ~Cf,o nE o»g 2 t/ s o* (q°, q, p)}^ 2 ) , 

(17) 

where _ 

Cdny = a fy—^ 1 (18) 

5 2/ 

with a given in Table [5] and f Y being a recoil factor [77j , which we approximate by 
/* ~ (1 - M D /M Y ). 

In Eqs. (|14p and f)lT[) . we include a form factor of monopole type at the D 
meson-baryon vertices by analogy to the one accompanying the Yukawa ttNN vertex 

IZZHSDI- 

A 2 

F L (q 2 ) = — ^ with A = 1.05 GeV. (19) 

A z + q z 

This form factor is suited for light hadrons, i.e., pion excitation of ph. However, 
it is unlikely that the range of A is the same when dealing with D mesons. There are 
indications that the form factor to account for off shell D mesons requires a value 
of A substantially larger [81]. We shall come back to this point at the end of the 
Results Section, reevaluating results with the heavy meson form factor and analyzing 
the uncertainties. 

With regards to the p-wave D + selfenergy, it turns out to be the same as for D° 
in symmetric nuclear matter p n = p p . 

For the D meson, we note that the p-wave D selfenergy would correspond to the 
diagrams in Fig. [3l which involve the difference between the sum of D and Y masses, 
and the nucleon mass. The contribution of those diagrams is negligible due to the 
large mass of the D and E c , A c . The same holds for the p-wave D selfenergy coming 
from the iV~ 1 S*(2520) excitation. 
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Figure 3: p-wave selfenergy diagrams of the D meson. 



3.3 The D meson spectral function 

The selfenergy of a I? meson in nuclear matter is given by the coherent sum of the 
s-wave and p-wave selfenergies: 

u D (q°,q, P ) =ngV,g>) + ngV,g>)+ngV,9>) ■ (20) 

Then, the D propagator is written in the medium as 

D D (q°,q,p) = 5 . (21) 

iq0) 2_^_ m 2 D _ UD{q0j ^ p) 

For later purposes, it is convenient to write the D propagator in terms of its Lehmann 
representation [82\ [83] 

Jo lq — uj + ir] q + w — ir\ ) 
where Sd and are the spectral functions of D and D, respectively, 

s D( D)(q ,q,p) = y - 1 2 • ( 23 

In the calculations we will ignore the selfenergy of the D. As in the case of the 
K with respect to the K, the D selfenergy is smaller than that of the D meson and, 
more importantly, it has no imaginary part from inelastic channels. Hence, it does 
not lead to modifications of the width of the states that we study here, which is the 
most striking change that we find. We will discuss this in more detail in what follows. 

P / 1 



/ * \ 

\ / 

\ P - q / 
D ^ — " 

Figure 4: (Color online) The DD loop function of the scalar meson. The shaded circle indicates 
the D selfenergy insertion. 
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q [GeV] 



q [GeV] 



Figure 5: (Color online) Real (left column) and Imaginary (right column) part of the s-wave 
D selfenergy for D three- momenta q = 0.15, 0.3 and 0.8 GeV as a function of the D energy q° at 
densities p = 0.5po and p = p$, with po = 0.17fm -3 the normal nuclear matter density. 




q [GeV] 



q [GeV] 



Figure 6: (Color online) Real (left column) and Imaginary (right column) part of the p-wave 
D selfenergy for D three- momenta q = 0.15, 0.3 and 0.8 GeV as a function of the D energy q° at 
densities p = 0.5/9q and p = p$, with po = 0.17fm~ 3 the normal nuclear matter density. 
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4 Two meson loop function in the medium 



The medium modifications are introduced in the two loop meson function by using 
the dressed two meson propagator in nuclear matter. As an example, let us evaluate 
the G function in the medium for the DD intermediate state: 



We shall dress the D propagator and leave the D propagator free. The reason to 
neglect the D selfenergy in the medium is that it is very small compared to its mass 
[56] . Indeed, the p-wave part is negligible as discussed at the end of Subsection 4.2. 
The s-wave part is equally small, but more importantly there is no absorption of D by 
nucleons, meaning that the DN does not decay to baryonic resonances which have c 
quarks instead of c . The analogy is clear with the K and K, where the K (analogous 
to D) can undergo absortion reactions KN —* 7rA, 7r£, while the K + (analogous to 
D) cannot be absorbed. Altogether this justifies to use the free propagator for D. 
We must evaluate the loop function of the diagram of Fig. [H where the blob in the 
D propagator symbolizes the D selfenergy insertion, indicating that we must use the 
D propagator in the medium. Then, we use Eq. (|22p for this propagator. Given the 
large mass of the D mesons, we can also neglect the part in the propagator of 
Eq. (|22p . which upon q° integration in the Eq. (|24p leads to a contribution of order 
of 1/(P° + 2co£)(q)), very small compared with the part coming from Sp. Thus, we 
can write 



Qmax = 0.85 and 0.9 GeV for the A(3700) and Z? s o(2317), respectively, equivalent to 
the use of dimensional regularization of [12] with the chosen substraction constants. 




Following the discussion in the former section, the DD loop function will appear 
in the case of the X(3700) hidden charm state. In the case of the -D s o(2317), we 
have K instead of D and, again, we use its free propagator neglecting the small 
structureless K selfenergy. 



In Figs. [5] and [6] we show the results of the s-wave and p-wave selfenergies of the 
D meson. We perform the calculations first by using the light meson version of the 
Dph form factor, Eq. (|19p . We will show results with a heavy meson form factor at 
the end of this section. In Fig. [5] the structures around energies 1.7 GeV and 2 GeV 
correspond to the excitation of the /iA c (2593) and /iS c (2800), where the A c (2593) and 
S c (2800) are l/2*~ dynamically generated resonances of the theory and h stands for 




(24) 





5 Results 



14 



hole of nucleon. With regard to the p-wave contribution, the structures seen in the 
p-wave selfenergy, with peaks for the imaginary part and oscillations in the real one 
around 1.4 to 1.8 GeV, correspond to the excitation of the different hY components, 
with Y = A c , E CJ S*. 

Note that the p-wave selfenergy is much smaller than the s-wave selfenergy, even 
up to momentum as large as 800 MeV/c. The main reason is that the s-wave selfen- 
ergy of the mesons goes roughly as the meson mass, from the Weinberg Tomozawa 
interaction, while the p-wave scales differently, roughly like the baryon mass. In the 
case of pions, the p-wave selfenergy was more important than the s-wave [M], but in 
the case of kaons the relative importance of the p-wave was already smaller |59[ 185] . 

In Figs. [7] and El we plot the G(P°,P,p) function for P = for the DD and 
DK loops for different densities. The effects of the density in the loop function are 
clearly visible in all cases. The imaginary part is largely increased at lower energies 
and collects strength below threshold of the DD or DK channels, respectively, due 
to the opening of new many body decay channels of the meson-meson system. As 
an example, let us take the DK loop function (see Fig. [9]). The D is renormalized 
and accounts for DN — * ttA c , 7rS c or DN — > A c , S c . Hence, the DK loop in the 
medium accounts for intermediate channels hirA c K, hir'E c K or hA c K, hTi c K which 
have a smaller mass than the DK system and open up at lower energies than DK 
threshold. The real parts of G are also sizeably modified around the thresholds as 
one can see in the figures. 
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Figure 7: (Color online) Loop function in the medium: Re G(P°, P, p) (left) and ImG(P°, P, 
(right) for DD, the channel with the largest coupling to the X(3700) meson. G{P ,P, p) is giv 
fromEq. (|25|) . 



In Fig. [10] we show the \T\ 2 for the D Q K + — > D°K + amplitude around the region 
of 2300 MeV for different densities. We can see that originally, at p = 0, the amplitude 
exhibits the pole of the 7J s o(2317), with zero width. This is the reason why it is out 
of the y-scale in the plot since \T\ goes up to infinite. As the density increases, 
we can see a slight shift of the mass, of the order of 15 MeV attraction at p = pq. 
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Figure 8: (Color online) Loop function in the medium: Re G(P°,P, p) (left) and ImG(P°, P, p) 
(right) for DK, the channel with the largest coupling to the D s o(2317) meson. G(P°, P, p) is 
given from Eq. (|25p . 



However, the increase in the width is more spectacular, which goes from zero in the 
free case to about 100 MeV at p = po, and 200 MeV at p = 2pQ. This is certainly 
a drastic relative effect, and even big in absolute value. The origin of the width in 
the medium is due to the opening of new channels DN — > A c , S c from the p-wave 
selfenergy and DN — > 7rA c ,7rS c from the s-wave selfenergy. On the other hand, the 
use of selfconsistency in the evaluation of the D selfenergy [551 ES] generates also some 
two nucleon induced D absorption channels like DNN — > NA C , irNA c , 7T-/VE C , etc. 
All these channels collaborate to make the D disappear inside the nuclear medium 
through DN or DNN inelastic reactions, where the D gets absorbed. 



D . _ 



N 



£ s0 (2317) 



s 
K 



Figure 9: Decay channel of the D s q(2317) in the nucleus into KttA c or KitYi c . 



The study of such decay channels in the nucleus would offer information on the 
coupling of the D s o(2317) to the DK, the basic building block of the resonance 
according to the underlying theory that we are using. As an example, the exploration 
of the decay channel of the -D s q(2317) in the nucleus into KirA c or Ktt'Ec channels, 
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Figure 10: (Color online) Z) s o(2317) resonance: \T\ 2 for the D°K + -> D°K+ amplitude for 
different densities. 

which corresponds to the cut in the diagram of Fig. [91 would provide combined 
information on the -D s o(2317) — ► DK coupling and the DN — > 7rA c (£ c ) cross section. 

The results with the large width for the .D s o(2317) are a consequence of the 
large coupling of this resonance to DK. This large coupling is guaranteed by the 
"compositness condition" of Weinberg [62j-[67j, as far as the resonance is dynamically 
generated and DK is the main building block. Should this resonance be a qq state 
as suggested in |86| or have any other structure, like qq with a mixture of DK 
components as suggested in [H71[S8], such large coupling would not appear [62] and, 
thus, the width in the medium would be much smaller than we predict here. Hence, 
investigating the widths in the medium provides extra information on the nature of 
these resonances. 

In Fig. [11] we show the same results for the hidden charm X(3700) resonance. 
The resonance begins with about 60 MeV in the real axis (36 MeV deduced from the 
imaginary part of the pole position) from its main decay into the 7777, 7777', and rfrf 
channels. The difference with [42], where the width was much smaller, is due to the 
fact that, here, we use the mixing of mesons ij and rf [60] . which makes the rj c a pure 
cc state rather than pure SU(4) state of the 15-plet of mesons. The study of this 
figure shows that at p = po the width has become as large as 250 MeV and, due to 
cusps on the multiparticle channels that open up, the strength of the resonance in the 
nuclear medium acquires a peculiar shape as the density increases. Here, the main 
decay channels in the nucleus are the same as in the case of the D s o(2317), replacing 
the K by a D. One should note that there is no relationship of the important decay 
channels to those that dominate the dynamical generation of the state. For example, 
the X(3700) couples mainly to DD but the DD threshold is above the X(3700) in 
free space and there is no decay into these components. The analogy can be made 
with the /o(980) which couples mostly to KK [891 EH] but decays into -kit in free 
space, since there is no phase space for KK decay. The coupling of the X(3700) 




17 



3.5x10 s 
3.0xl0 5 
2.5x10 s 
2.0xlO S 
1.5xl0 S 
1.0x10 s 
S.OxlO 4 



3200 3400 3600 3800 4000 

s" 2 [MeV] 

Figure 11: (Color online) X(3700) resonance: \T\ 2 for the D°D° -> D°D° amplitude for 
different densities. 

to the 7777, 777/ etc components is small, reflecting the fact that these channels are 
relatively unimportant in the structure of the resonance. However, because there 
is large phase space for these decays, they are mostly responsible for the X(3700) 
width. 

In the medium things might be different. Indeed, in the presence of nucleons 
we can have DN — > 7rA c , 7rX c as we mentioned above, and new decay channels for 
the X(3700) appear, such as X(3700)iV DttA c , D-kY> c [see Fig. El (changing K 
by D)], which have nearly 400 MeV phase space available. The strong coupling 
of the X(3700) to DD and the opening of these new decay modes makes now the 
width in the medium sizeable. We note that the trend of changes of \T\ 2 with nuclear 
density is typically observed in resonance properties with temperature and/or density 

[an si [93]. 

As mentioned at the end of Subsection 4.2, in order to analyze the uncertainties 
linked to the use of different form factors, we reevaluate the results using a form 
factor for the off shell D mesons which we obtain from the work of [81]. The heavy 
meson form factor is now 

A 2 2 

F H (q 2 ) = R ~ - with A D = 3.5 GeV. (26) 
A Z D - q 2 

The value of q 2 in Ffj(q 2 ) is taken for the configurations that give rise to the width 
in diagram of Fig. [9] when the D selfenergy comes from particle-hole (ph) excitation 
of Fig. [2j This means that, for the case of the D s q(2317), we place the K and A c h, 
T, c h or E*/i on shell. This occurs for any value of the running variable q at a value 
of q° given by 

„o w . . / „ \ - fh „ A 1 / 2 ((M Ds0 +m jV ) 2 ,M 2 ,m^) 

q - M Ds0 - u K {q on ) with q on - 2 (M Da0 +m N ) (2?) 
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p 


M[MeV] 


T[MeV] 


0.0 


2316±5 





0.5 


2306±17 


58±10 


1.0 


2295±23 


115±25 


1.5 


2283±25 


150±25 


2.0 


2274±31 


190±30 



Table 6: Mass and width for the _D s0 (2317) at different densities with error bands due to 
the uncertainties of our model. 



p 


M[MeV] 


T[MeV] 


0.0 


3710±18 


60± 10 


0.5 


3691± 10 


135±20 


1.0 


3638± 15 


255±25 


1.5 


3599±15 


320±25 


2.0 


3565±29 


340±25 



Table 7: Mass and width for the X(3700) at different densities with error bands due to 
the uncertainties of our model. 

where M can be M\ c , My, c or depending on the type of ph excitation. We 

anticipate that, given the minor relevance of the p-wave selfenergy, the differences 
with respect to the former calculation using the light meson form factor will be small. 
This is indeed the case as can be seen in Fig. [T2| where we show |T| 2 for the case 
of the D s q(2317) for the two types of form factor. As we can see there are small 
differences in the position of the peak and no difference in the width. 

As indicated before, in Figs. [10] and [11] we have shown the squared amplitudes 
for different densities to facilitate the comparison. As indicated at the end of Section 
2, we also evaluate the uncertainties in the results. We show in Tables [6] and [7] the 
results obtained for the mass and width of the resonances at different densities, with 
their uncertainties obtained from a Monte Carlo sampling of the parameters f n and 
fjj as indicated in Section 2 and evaluated from the plots of |T| 2 in the real axis. 
The results in these tables are all evaluated using the heavy meson form factor of Eq. 
(|26p . The relevant information from these tables is that the differences between the 
widths at p = or p = po for both resonances are much bigger than the uncertainties 
in the width from uncertainties in the model. In the case of the masses, we do not 
see appreciable shift in the mass compared to the uncertainties for the case of the 
£> s0 (2317), see Tabled For the X(3700), the shift from p = to p = p is of the 
order of 70 MeV, see Table [7j smaller than the width in the medium. Hence we 
cannot make any strong point concerning a possible shift of the masses. 

In view of future experiments measuring medium modifications of these reso- 
nances, we can recall the method that has proved most efficient in measuring nuclear 
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Figure 12: (Color online) Comparison of \T\ 2 in the case of the D s0 (2317) resonance for 
the two different form factors at p — p : type 1 (dashed line) with F H (q 2 ) of Eq. (T2T)j) and 
type 2 (solid line) with F L (q 2 ) of Eq. <HM . 



widths: the transparency ratio. The direct measurement in experiments of the in 
medium increased width is not easy because the observed decay channels usually 
come from the resonances that have escaped from the nucleus, so the density at the 
decay place is zero or very small [2j. Therefore, one should look at the production rate 
as a function of the mass number normalized to a particular nucleus (transparency 
ratio). This magnitude, which measures the survival probability, is very sensitive to 
the absorption rate of the resonance inside the nucleus, i.e., the in medium resonance 
width. This procedure has been succesfully used for the <f> and u> production in nuclei 
in [U [94] with the help of relatively easy tools of analysis j3j [95] . 



6 Conclusions 

We have evaluated the selfenergy of low lying scalar mesons with open and hidden 
charm in a nuclear medium, concretely of the -D s o(2317) and the theoretical hidden 
charm state X(3700). The many body calculation has been done following the lines 
of previous studies in the renormalization of the light scalar mesons in the nuclear 
medium. The medium effects for the _D s o(2317) and X(3700) resonances are spectac- 
ular. Those resonances, which have zero and small width in free space, respectively, 
develop widths of the order of 100 and 200 MeV at normal nuclear matter density, 
respectively. The study also allowed us to trace back the reactions in the medium, 
which are responsible for the decay width of these mesons and which could be inves- 
tigated in future reactions at hadron facilities. The option of looking at transparency 
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ratios was also suggested as a mean to investigate the widths of these mesons in nu- 
clei. It was also discussed that the experimental study of this width and the medium 
reactions contributing to it provide information on the basic features of the resonance 
and the selfenergy of the D meson in a nuclear medium. In other words, the experi- 
mental analysis of those properties is a valuable test of the dynamics of the D meson 
interaction with nucleons and nuclei, and the nature of the charm and hidden charm 
scalar resonances, all of them topics which are subject of much debate at present. 
The results obtained here should stimulate experimental work in hadron facilities, 
in particular at FAIR [96] . where the investigation of charm physics is one of the 
priorities. 
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